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ABSTRACT 

We present a study of the circumnuclear region of the nearby Seyfert galaxy Mrk 573 using Chandra, 
XMM-Newton and Hubble Space Telescope (HST) data. We have studied the morphology of the soft 
(< 2 keV) X-rays comparing it with the [O III] and Ha HST images. The soft X-ray emission is 
resolved into a complex extended region. The X-ray morphology shows a biconical region extending 
up to 12 arcsecs (4 kpc) in projection from the nucleus. A strong correlation between the X-rays and 
the highly ionized gas seen in the [O III]A5007 A image is reported. Moreover, we have studied the line 
intensities detected with the XMM-Newton Reflection Grating Spectrometer (RGS) and used them 
to fit the low resolution EPIC / XMM-Newton and ACIS/ Chandra spectra. The RGS/ XMM-Newton 
spectrum is dominated by emission lines of C VI, O VII, O VIII, Fe XVII, and Ne IX, among others 
highly ionized species. A good fit is obtained using these emission lines found in the RGS/ XMM- 
Newton spectrum as a template for Chandra spectra of the nucleus and extended emission, coincident 
with the cone- like structures seen in the [O III] /Ha map. The photoionization model Cloudy provides 
a reasonable fit for both the nuclear region and the cone-like structures showing that the dominant 
excitation mechanism is photoionization. For the nucleus the emission is modelled using two phases: 
a high ionization [log(U)=1.23] and a low ionization [log(U)=0.13]. For the high ionization phase the 
transmitted and reflected component are in a ratio 1:2, whereas for the low ionization the reflected 
component dominates. For the extended emission, we successfully reproduced the emission with two 
phases. The first phase shows a higher ionization parameter for the NW (log(U)=0.9) than for the 
SE cone (log(U)=0.3). Moreover, this phase is transmission dominated for the SE cone and reflection 
dominated for the NW cone. The second phase shows a low ionization parameter (log(U)=-3) and is 
rather uniform for NW and SEconcs and equally distributed in reflection and transmission components. 
In addition, we have also derived the optical/infrared spectral energy distribution (SED) of the nucleus 
from high spatial resolution images of Mrk 573. The nuclear optical/infrared SED of the nucleus has 
been modeled by a clumpy torus model. The torus bolometric luminosity agrees very well with the 
AGN luminosity inferred from the observed hard X-ray spectrum. The optical depth along the line of 
sight expected from the torus modeling indicates a high neutral hydrogen column density in agreement 
with the classification of the nucleus of Mrk 573 as a Compton-thick AGN. 

Subject headings: galaxies:active - galaxies: nuclei - galaxies :Scyfcrt - galaxies: individual (Mrk 573) - 
infrared : galaxies 



1. INTRODUCTION 

Visible signatures of the direct interaction be- 
tween the active galactic nuclei (AGN) and their 
host galaxies include kpc-scale [0 III]- and Ha- 
cmitting regions, the so-called extended narrow-line 
region (ENLR), which is observed in many nearby 
Seyfe rt galaxies dSchmitt et al.l 120031 iVeilleux et al 



20031 iWhittle fe Wilson! 12001 iRamos Almeida et al 



2006( 1. Understanding the AGN-host galaxy interaction 



and feedback is cr ucial for the study of both galaxy 
and AGN evolution (ISilk fe Reeslll998l:lKauffmann et all 
[20031: IBopkins et al.ll2006t ISchawinski et al. 20071 120091k 

One of the most promising ways to study these in- 
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terac t ions is through soft X-rays (|Crenshaw et al.|["l999l . 
[20031: iDai et al J 120081) . In the unified picture of AGN 
(|Antonuccilll993l ) the soft X-ray spectra of Type-2 Seyfert 
galaxies are expected to be affected by emission and scat- 
tering from the medium, which is also strongly influenced 
by the nuclear continuum. However, soft X-rays can also 
be produced through mechanical heating, as in shocks 
driven by supernova explosions in nuclear star-forming 
regions. Indeed, it is quit e plausible that both effects 
are important (e.g. Mrk 3: iSako et al.ll200fl ). This soft 
X-ray emission provides the opportunity to obtain an 
X-ray diagnostic for the physical pr operties of the in- 
teracting interste l lar medium [ISM ] (lYoung et al.l 120011: 
l Yang et all [2001: jO^ielFall 120001 120031: iBianchi et al.l 
120061 : lEvans et al J 120011 iKraemer et al.ll2008l ). 

The galaxy Mrk 573 has been extensively studied 
by many authors. Its active nucleus is hosted in 
an (R)SAB(rs)0+ galaxyQ Mrk 573 is well-known for 
its extended, richly structured circumn uclear emission- 
line regions (jTsvetanov fe Walshl [l£92). It has long 

6 T he redshift of the galaxy has been measured to be z = 
0.017 (Ruiz et al. 2005), which implies a physical scale of 333 
pc arcsec -1 , using Hq = 75 km s _1 Mpc - 1 . 
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been known to show a biconical st ructure and bright 
arcs and knots of l ine-emitting gas (|Ferruit et all 119991 : 
iQuillen et al.l I1999D that are strongly aligned and in- 
teracti ng with a kiloparsec- s cale l o w-power radio out- 
flow dPogge k De Robertisl 119931: iFalcke et al.1 119981 : 
iFerruit et all 119991 ) . iSchlesinger et al.l (|2009( ) have re- 
cently studied the STIS/ HST spectrum, finding that the 
ENLR optical spectrum is consistent with photoioniza- 
tion by the AGN. 

In their study of the A CTS / C handra spectrum 
of Mrk573 iGuainazzi et al.l (|2005l ) showed that it 
is consistent with a Compton-thick source (i.e. 
Nh > 1.6 x 10 _24 cm~ 2 ), showing a large EW of the 
FeKa emission line and a steep photon spectral in- 
dex (r = 2.7±0.4). These characteristic X-ray prop- 
erties, together with high-quality LIRIS near-infrared 
spectroscopy, allowed this galax y to be reclassified as 
an obscured narrow-line Seyfert 1 (|Ramos Almeida et al.l 
l20Ml2009ah . 

RGS / XMM-Newton high resolution observations of 
M rk573 were reported a mong a sample of 69 objects 
bv lGuainazzi et al~l (|2008|) . whose spectrum shows strong 
emission coming from the O VII triplet and Lya O VIII 
features (see their figure 1), but only few emission line 
fluxes were reported in their work. However, the anal- 
ysis of the extended soft X-ray emission had not been 
previously reported in the literature. 

In this paper, we discuss the extended kpc-scale emis- 
sion of Mrk 573 using ACTS / Chandra high resolution im- 
ages, RGS / XMM-Newton high resolution spectra and 
HST imaging. While RGS/ XMM-Newton data give the 
opportunity to use emission line diagnostics to under- 
stand the nature of this emission, Chandra high reso- 
lution images give the opportunity to spatially resolve 
the emission and study their properties separately. HST 
images allow us to establish the connection between the 
optical structure and the soft X-ray emission. The pa- 
per is presented as follows. Section 2 describes the data 
reduction and processing. We present the study of the 
circumnuclear morphology and the X-ray spectral anal- 
ysis in Sections 3 and 4, respectively. The origin of this 
soft X-ray emission is discussed in Section 5 and the op- 
tical/infrared spectral energy distribution is studied in 
Section 6. Finally, conclusions and the overall picture 
for Mrk 573 are presented in Section 7. 

2. OBSERVATIONS AND DATA PROCESSING 
2.1. XMM-Newton data 

We retrieved from the HEASARC archive the XMM- 
Newton observation of Mrk 573 taken on 2004 January 15 
(ObsID 0200430701). RGS data were processed with the 
standard RGS pipeline processi ng chains incorpora ted in 
the XMM-Newton SAS v.8.0.1 (| Gabriel et alll200l . Dis- 
persed source and background spectra (using blank field 
event lists) were extracted with automatic RGS extrac- 
tion tasks. The net exposure time is 9 ks after flare 
removal and the net count rate is 0.21 count s 

In this paper we also us ed lower resolution s pectrum 
from the EPIC pn camera (jStruder et al.ll2001l) . Source 
regions were extracted within a circular region of 25 arc- 
seefl radius centered at the position given by NEE0. We 

7 This radius contains the 80% (85%) of the Point Spread Func- 



also used eregionanalise SAS task to compute the best 
centroid for our source. The difference between the best 
fit centroid and the NED position is 2 arcsec (P.A. 126°). 
It implies ~1% of error on the final flux, according to the 
PSF of EPIC pn instrument. Background was selected 
from a circular region in the same chip as for the source 
region, excluding point sources. Regions were extracted 
by using the evselect task and pn redistribution ma- 
trix, and effective areas were calculated with the RMF- 
GEN and ARFGEN tasks, respectively We also binned 
the EPIC/pn spectrum to give a minimum of 20 counts 
per bin before background subtraction to be able to use 
the x 2 as the fit statistics using the grppha task. Note 
that the background is only 2.4% of the total number of 
counts (~ 50 counts). Thus, the spectral binning can be 
done before background substraction in this source. 

2.2. Chandra data 

Mrk 573 was observed by Chandra on 2006 Novem- 
ber 11. Level 2 event data from the ACTS instru- 
ment were extracted from the Chandra archive^ (ObsID 
7745). The data were reduced with the CIAO 3.43 data 
analysis system and the Chandra Calibration Database 
(caldb 3. 4.07^1). The exposure time was processed to ex- 
clude background flares using the LCLCLEAN.SL taslQ hi 
source-free sky regions of the same observation. The net 
exposure time after flare removal is 35 ks and the net 
count rate is 1.5 s . The nucleus has not significantly 
piled up. 

Chandra data include information about the photon 
energies and positions that was used to obtain energy- 
filtered images and to carry out sub-pixel resolution spa- 
tial analysis. Although the default pixel size of the Chan- 
dra/ ACTS detector is 0.492 arcsec, smaller spatial scales 
are accessible as the image moves across the detector pix- 
els during the telescope dither, therefore sampling pixel 
scales smaller than the default pixel of Chandra/ ACIS 
detector. This allows sub-pixel binning of the images. 
Similar techniques were applied for the analysis of Chan- 
dra observations of, for example, the SN1987A remnant 
([Burrows et al.ll2000D . 

In addition to the high-spatial resolution analysis, we 
applied smoothing techniques to detect the low-contrast 
diffuse emission. We applied the adaptive smoothing 
CIAO t ool CSMOOTH based on the algorithm devel- 
oped bv lEbeling et all (|2006l ). CSMOOTH is an adaptive 
smoothing tool for images containing multiscale complex 
structures and preserves the spatial signatures and the 
associated counts, as well as significance estimates. A 
minimum and maximum significance S/N level of 3 and 
4, and a scale maximum of 2 pixels were used. 

We also performed spectral analysis of the nuclear and 
extended emission using CIAO software. Background re- 
gions were defined by source-free apertures around the 
soft X-ray emission regions. Response and ancillary re- 
sponse files were created using the CIAO mkacisrmf 

tion (PSF) at 1.5 keV (9.0 keV) for an on-axis source with EPIC 
pn instrument. 

8 http:/ /nedwww. ipac.caltech.edu/ 

9 http://cda.harvard.edu/chaser/ 

10 http://asc.harvard.edu/ciao 

11 http://cxc.harvard.edu/caldb/ 

12 http: / /cxc. harvard.edu/ciao/download/scripts/ 
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TABLE 1 

High spatial resolution nuclear fluxes 



1000 



A .ill 1 vi 1 
A re f \I- Lln ) 


F (m lv) 


Tn <3l~vi imont" /hiI+pi* 


Rcf. 


0.675 


0.002 


WFPC2/F675W 


a 


0.814 


0.003 


WFPC2/F814W 


a 


1.100 


0.036 


NICMOS/F110W 


a 


1.650 


0.479 


NICMOS/F160W 


a 


2.120 


3.20 


NSFCam/K' 


b 


3.510 


18.8 


NSFCam/L 


b 


4.800 


41.3 


NSFCam/M 


b 


10.36 


177 


T-ReCS/N 


c 


18.30 


415 


T-RcCS/Qa 


c 



Note. — Errors in flux densities are ~10% for HST fluxes, 
~20% for the 3 m NASA IRTF, ~15% for T-ReCS/N. and -25% 
for T-ReCS/Qa. 

References. — (a) this work ; (bl lAlonso-Herrero et all H2003T ); 
fcl lRamos Almeida etaTI <2009bl ). 



and MKWARF tools. To be able to use the x 2 as the fit 
statistics, the spectra were binned to give a minimum of 
20 counts per bin before background subtraction using 
the GRPPHA task, included in FTOOLS0 
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2.3. Optical and near-IR data 

In order to perform a detailed comparison between X- 
ray and optical cone-like structures we have retrieved 
HST /WFPC2 narrow-band optical images, centered at 
5343 A([0 III]) and 6510 A(Ha). These images were ob- 
tained within the Cycle 4 program 6332 (PI Wilson) and 
retrieved from the Hubble Leg acy Archive. Details of the 
observations can be found in iFalcke et all (|1998f ). The 
two images were re-centered assuming that the peak in 
both images corresponds to the same locus (|Ferruit et all 
119991 ). The same assumption was used to align opti- 
cal and X-rays images. Note that high accuracy is not 
needed due to the lower Chandra resolution. 

We have also investigated the nuclear properties of 
Mrk 573 in the optical and near-IR ranges, based on high 
spatial resolution broad-band images. We have found 
optical and near-IR broad-band images taken with the 
WFPC2 and NICMOS cameras on board the HST. In 
particular, we used F675W and F814W WFPC2 images 
(program 5746; PI Machetto), and F110W and F160W 
NICMOS images (program 7867, PI Pogge). All data 
were retrieved from the Hubble Legacy Archive. Post- 
pipeline images have been cl eaned of cosmic ra ys us- 
ing the IRAF task lacos_im (jvan Dokkuml l200lT ) . For 
the analysis we have first separated the nuclear emission 
from the underlying host galaxy emission. We have ap- 
plied the two-dimensional i mage decomposition GALFIT 
program (|Peng et al.ll2002l ) to fit and subtract the unre- 
solved component (PSF) . Model PSFs were created using 
the Tiny Tim software package, which includes the optics 
of HST plus the specifics of the camera and filter system 
(|Kristlll993T ) . The filters used here are not contaminated 
by line emission except for the case of F110W, which 
contains strong Pa/3 emission. We have used the nuclear 
near-IR spectrum of Mrk 573 obtained with th e LIRIS 
spectrograph (see IRamos Almeida et al.l I2009a[) to sub- 
tract this emission line contribution. The Pa/3 nuclear 
emission (including the broad and narrow components) 

13 http:/ /heasarc. gsfc.nasa.gov/docs/software/ftools 



Fig. 1. — Radial profile of the soft (0.2-2.0 keV) X-ray emission 
of Mrk 573 (black dotted). Continuous linos correspond to the fit 
to the inner (< 1.5arcsecs) and outer parts of the radial profile. 
Dashed line shows the radial profile of the simulated PSF at ~1 
keV. The model used to fit the radial profile is of the form: y(x) = 

A X e~ x / 2B . An additive constant was added to the outer parts 
of the radial profile fit to include the minimum S/N level (dotted 
line). The FWHM of the inner parts is 0.59 arcsecs, consistent 
with the FWHM of the PSF (0.43 arcsecs). 

amounts 3.17 x 10~ 14 erg cm~ 2 s _1 . The corrected values 
of the flux obtained for the nuclear component in each 
filter are given in Table [T] 

3. CIRCUMNUCLEAR MORPHOLOGY 

Two broad-band images using Chandra/ ACIS dat;Q 
were created: (1) soft band at 0.2-2 keV, and (2) hard 
band at 2-10 keV. PSF simulations were carried out us- 
ing information on the spectral distribution and off-axis 
location of the system as inputs to ChaRT PSF simula- 
tor[3 Hard band (>2 keV) shows a point-like morphol- 
ogy, consistent with the FWHM PSF simulations, while 
soft band (<2 keV) shows an extended morphology (see 
Figure [J) . 

The soft X-ray image shows a complex extended emis- 
sion with a bipolar structure aligned along PA~ 122° 
(Figure It extends 8.5 arcsec to the NW and 6.5 
arcsec to the SE. As already mentioned, Mrk 573 shows 
a biconical morphology in the [O III] emission, which 
resembles that observed in the soft X-rays. 

The ratio of [OIII] to H recombination lines is com- 
monly assumed to be an indicator of the ionization de- 
gree in line emitting photoionized gas. We have com- 
puted the [O III] /Ha map by dividing the flux-calibrated 
HST/WFPC2 images in those pixels where the signal-to- 
noise is higher than 3. A smoothing algorithm using the 
median within a box of 0.5 arcsec was then applied. Note 

14 We used Chandra data for morphological analysis since it 
provides the best X-ray spatial resolution. 

15 http:/ /asc. harvard.edu/chart/index. html 
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that this ratio image is only used for morphological com- 
parison and not for any quantitative analysis since the 
Ha image is contaminated with [N II] emission. The use 
of Ha instead of H/3 to compute the ionization map could 
be affected by reddening variations, as shown close to the 
nucleus of Mrk 573. However, the extinction map (Fig. 
[3]) shows a very different morphology to that observed in 
the ratio [O III] /Ha map. This implies that the basic 
morphology of the optical ionization map cannot be ex- 
plained by differential extinction, despite a certain frac- 
tion of the emission line map variations could be due to 
it. 

The X-ray arc-like structure at 10 arcsec to the SE 
resembles the [O III] emission. However, the X-ray emis- 
sion does not show the bridge between this structure 
and the inner parts of the SE cone seen in the [O III] 
image (also seen in [O III] /Ha). The NW structure is 




-4 -2 2 4 

RA (arcsec) 



Fig. 3. — I-H color map of Mrk573. Left is the East and up is 
the North coordinates. 



coincident with the [O III] emission and strongly simi- 
lar to the [O III] /Ha morphology (see Figure [2]). These 
facts indicate a link between the soft X-ray emission and 
the optical ionized gas, although the detailed structure 
would depend on the small-scale ionization structure of 
the medium. 

We have also explored the circumnuclear extinction. 
For this purpose we have used the broad-band images 
after subtracting the unresolved component to build a 
colour map of the circumnuclear region of Mrk 573. The 
colour map was created using the F814W and the F160W 
images, which are equivalent to the / and i?-bands, re- 
spectively. These filters were selected because they are 
not critically contaminated by emission lines , contr ary to 
the colour maps presented bv lQuillen et al.l (|1999f ) . The 
I — H map covering the central 10 arcsec is shown in 
Figure [3l The reddest (or darkest in black/white) re- 
gions show the location of the strongest obscuration, or 
equivalently, the highest dust concentration. The colour 
map reveals a dust lane crossing the nucleus in the N- 
S direction, which bends after 2" resembling incipient 
spiral arms. The axis delineated by this structure is ori- 
ented at about 5 5° with respect to t he alignment of the 
radio structure (jFalcke et al.1 119981 iKinnev et all 120001 
PA ~ 125°). The radio and the dust lanes axis should 
be almost perpendicular if the dust lanes were the outer 
parts of the postulated dusty torus for Type-2 nuclei. 
However, this is not necessarily the case, since the dust 
lane is observed at much larger scale than that of the 
postulated torus, changes in the dust plane may take 
place when approaching the inner nucleus. On the other 
hand, projection effects could mas k the actual orienta- 
tion o f both structures, for instance iTsvetanov fc Walsh! 
(|1992f l proposed that the ionization cone is very inclined 
(35°) with respect to the plane of the sky. In addition, 
there is also a tongue extending 0'.'5 towards the N-NW 
which ends in two small blobs that resembles a struc- 
ture observed in the excitation maps (see Figure [2]) • The 
presence of dust within the ioniz ation cone has be en re- 
ported before only in NGC 1068 (|Bock et al.ll200Q[ ). To- 
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TABLE 2 

Emission lines detected in the RGS / XMM-Newton spectrum. 



Name 




Energy 


Flux 




(A) 


(keV) 


10" 5 



CV He 7 


33 


,426 





,376 


10 7 27 - 3 
±u - ' 0.0(1 


NVI Hea (f) 


29 


,534 





,420 





on 827 
uu o.o6 


NVI Hea (i) 


29 


,083 





,426 


2 


64 57.6 

on 185 

uu 0.00 


NVI Hea (r) 


28 


,787 





,431 





CVI Ly/3 


28 


,466 





,436 


3 


18 7.10 

HI 


NVII Lya 


21 


,781 





,500 


1 


03 2 51 


OVII Hea (f) 


22 


101 





,561 


5 


OVII Hea (i) 


21 


,803 





,569 





■oog : g§ 


OVII Hea (r) 


21 


,602 





,574 


3 


76 7 - 58 


OVIII Lya 


18 


,969 





,654 


2 


Ll 1.00 


OVII He7 


17 


,768 





,698 


1 


60 Hf 


FeXVII 3s-2p 


17 


,078 





,726 





7"> 1:84 
,o 0.00 


OVII RRC 


16 


,771 





,739 


1 


07 2 - 65 
u ' 0.03 


FeXVII 3d-2p 


15 


,010 





,826 





fit; 1.45 


OVIII RRC 


14 


,228 





,882 





4gl'§3 


NeIX(f) 


13 


,698 





,905 





90 8:82 

i2?:?8 
lz 0.00 


NeIX(i) 


13 


,552 





,915 





NeIX(r) 


13 


,447 





,922 





76 1 - a7 


a Units 


arc 


ph cm 


2 s" 1 . 







wards the nucleus of the galaxy there is an excess of 
reddening which can be attributed to a natural increase 
in the extinction due to higher dust concentration. As- 
suming an intrinsic colour sim ilar to that observed in the 
disk of galaxies (I - H - 2; IMoriondo et all 119981 ). we 
have estimated a value of Ay = 6.5 mag, which results 
in N H = 1.2 x 10 22 cm~ 2 . 



4. X-RAY SPECTRAL ANALYSIS 



X-ray spectral analysis of the observed soft X-ray ex- 
tended emission is crucial to determine the excitation 
mechanism of the plasma, and its relationship to the op- 
tical bicone-like structure (see previous section). The 
combination of RGS / XMM-Newton high spectral resolu- 
tion and ACIS/ Chandra high spatial resolution data is 
key to achieve this purpose. In this section we describe 
in detail the methodology and main results obtained. In 
Section [5] we discuss the origin of this extended emission 
based on the results presented in this section. The anal- 
ysis of the spectral counts was performe d using the sof t- 
ware package XSPEC (version 12.4.C0; IArnaudl[l996]) . 

4.1. RGS '/XMM-Newton high resolution spectra 

Soft X-ray emission in Seyfert galaxies has been proven 
to consist of a plethora of emission lines plus a small frac- 
tion of continuum em ission that can be des cribed with a 
single flat power-law (|Guainazzi et al.ll2008[ ) with a fixed 
spectral index of T = 1. We obtained the emission line 
fluxes of the central 30 arcsec region (note that this in- 
cludes the nucleus and circumnuclear emission) using the 
RGS / XMM-Newton data. We searched for the presence 
of 37 emission lines of C, O, N, Si, Mg and Fe species 
by fitting the spectra of the two RGS cameras to Gaus- 
sian profiles together with a continuum. We used Cash 
statistic for this purposes. 

The triplet fits were performed keeping the relative dis- 
tance between centroids in energy and the centroid en- 

16 http:/ /cxc. heasarc.gsfc.nasa.gov/docs/xanadu/xspec/, 



ergy was left as a free parameter. A line was considered 
detected when the flux was higher than at the la level. 
The resulting RGS spectrum and detected emission lines 
are presented in Figured] and Table [2] respectively. All 
energy centroids are c onsistent with th e labor atory value 
given the error bars. IGuainazzi et al.l (|2008|) previously 
studied the RGS/ XMM-Newton spectra of Mrk 573. Un- 
fortunately, they only reported some of the lines, all of 
them agreeing with our emission line fluxes. The most in- 
tense emission lines comprising the RGS / XMM-Newton 
spectrum are: C VI Ly/3, O VII (r), O VII (f), O VIII 
Lya. O VII H7, O VII RRC, Fc XVII 3d-2p, and Ne IX 
(r). ' 

4.2. EPIC '/XMM-Newton low resolution spectra 
The fit of EPIC/ XMM-Newton data with a thermal 
model produces poor results below 2 keV (\ 2 ~16). In- 
stead, a model composed of multiple emission lines was 
tried. Taking advantage of the RGS/ XMM-Newton fit, 
we imposed that the intensity of the lines in the low res- 
olution spectra fit do not exceed the RGS measurements. 
This is acceptable because the cross-calibrations between 
EPIC and RGS instruments shows a no rmalization con- 
stant in the range of 0.9 to 1.0 (see IPlucinskv et al.l 
120081 ). The assumed Gaussian width is 100 eV. Note 
that for EPIC (and also ACIS / Chandra) data we do not 
question the existence of the emission lines detected on 
the RGS / XMM-Newton but we use them as a template. 
Triplets were fitted using the total flux of all components 
of the He-like lines O VII, N VI, and Nc IX. The contin- 
uum emission was fitted to a power-law to be consistent 
with the high spectral resolution analysis. However, this 
fit has poor statistics (x 2 > 2). Five lines were added at 
energies above 0.95 keV in order to achieve an acceptable 
fit(x 2 = 0.8): FeXXat0.97keV(x 2 =5.64),NeXLyaat 
1.02 keV (x 2 =5.62), NcIX He S at 1.16 keV (x 2 =3.72), 
MgXI triplet at -1.33 keV (x 2 =0.82), and Si XIII 
triplet at 1.84 kcV(x 2 =0.80). The final fit is shown 
in Figured The low spectral resolution EPIC / XMM- 
Newton spectrum shows the following intense emission 
lines: CV HC7, CVI Ly/3, NVII Lya, OVII triplet, 
OVIII Lya, O VII He7, OVII RRC, FeXVII 3d-2p, and 
NelX triplet, NeX Lya, NelX He <5, and Mg XI triplet. 
In the best-fit model the flux of the OVIII RRC feature 
appears negligible and the adjacent line FeXVII 3d2p 
is present. This is in contrast to what happens in the 
ACIS /Chandra nuclear spectrum (see Section 14.31 and 
Table [3]). In order to check the compatibility of the re- 
sults of both spectra, we have imposed a zero intensity 
to the line FeXVII 3d2p finding a good fit. In this case, 
the OVIII RRC feature shows a similar flux compared 
to the flux reported using the Chandra nuclear spec- 
trum. It is very likely that both features share the flux as 
measured in the RGS / XMM-Newton spectrum, although 
they cannot be distinguished in the low resolution spec- 
tra. The best-fit model has an absorbed 0.5-2.0 keV flux 
of Fo. 5 _2.o kcv =3.38 (3.37-3.83) xl0~ 13 erg s" 1 cm- 2 , 
corresponding to an unabsorbed rest frame luminosity 
of Lo.5-2.0 kov =2.0 (1.9-2.2) xlO 41 erg s" 1 . We have 
only included Galactic absorption in our model, which 
corresponds to Nh = 2.52 x 10 20 . 

4.3. ACIS /Chandra spectra 
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Fig. 4. — High resolution spectrum of Mrk573 obtained with RGS/ XMM- Newton instrument. 
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Fig. 5. — Spectral fits (top panel) and residuals (bottom panel) 
for the nuclear spectrum of Mrk 573 using low resolution EPIC 
pn/ XMM-Newton data. Dotted-lincs show the Gaussian compo- 
nents used to fit the spectrum. 



We cannot separate the contribution of NW and SE re- 
gions to the RGS or EPIC/ XMM-Newton spectra. How- 
ever, this is possible by using the lower spectral resolution 
but better spatial resolution of ACIS / Chandra data. 

We extracted spectra from the nucleus (R = 1 arcsec) 
and two conical regions coincident with the extension of 
the emission, as seen in Figure [51 All extraction regions 
were centered at (RA, Dec)=01:43:57.78,+02:20:59.4. 
For the conical regions, we used an annulus, centered 



at the same position than the nucleus, with an inner 
and outer radius of 1.5 and 5 arcsec, respectively. The 
cones are defined by an opening angle of 60° centered at 
PA = 325° (cone NW) and at PA = 145° (cone SE). 

The study of the emission above 2 keV is beyond the 
scope of this paper since the nucleus do minates there and 
it ha s been already well studied before (jGuainazzi et al.l 
120051 ) . Therefore, channels above ~2 keV were ignored 
in the spectral fit. Again, the thermal model gives a poor 
fit with some residuals below 2 keV (\ 2 ^3 for both cone 
regions) . 

We used the same model reported in the RGS / XMM- 
Newton spectrum assuming the emission line fluxes found 
in that case, as upper limits to fit the nucleus, and 
the NW and SE cone-like structures (the same as the 
EPIC / XMM-Newton data mentioned in the previous 
section). Moreover, the power-law component was re- 
moved in the cone-like structures because we expect this 
component to be detectable only in the nuclear region. 
The assumed Gaussian width is 100 eV. The spectra are 
shown in Figure [7] and final emission line fluxes are re- 
ported in Table [3] 

Most of the emission lines in the EPIC / XMM-Newton 
spectrum are also detected in the nuclear spectrum us- 
ing Chandra data (CV FPy, N VII Lya, O VII triplet, O 
VIII Lya, Ne IX triplet, Ne X Lya, and Mg XI triplet). 
Five of them were present only in the EPIC/ XMM- 
Newton spectrum (C VI Lya, O VII FPy, O VII RRC, 
Fe XVII 3d-2p, and Ne IX He<5). We note that the flux 
measurements in the ACIS/Chandra spectrum are com- 
patible with those of the EPIC/XMM- Newton includ- 
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Fig. 6. — Soft (0.2-2 keV) X-ray Chandra image with the two 
extracted regions (NW and SE) overplotted. 



ing cases where only upper limits can be estimated in 
one of the spectra. In contrast with the EPIC/ XMM- 
Newton spectrum, the Ne IX HeS line at 1.16 keV was not 
needed, instead two more lines were added (NelX HC7 
and NeX Ly-7 at 1.13 and 1.22 keV, respectively) in or- 
der to get the best fit (x 2 =l-l). Note that in the case 
of Ne X Ly7 there are several transitions of Fe XX that 
could be contributing to this line. The inclusion of these 
two lines (i.e. NelX HC7 and NeX Ly7) instead of 
the NelX KeS line also provides an aceptablc fit to the 
EPIC / XMM-Newton spectrum. The line intensities are 
reported in Table [3] (within brackets). The final fit is 
equaly acceptable. A similar case is affecting the OVIII 
RRC and the FeXVII 3d-2p lines (see Section OJ . 

The emission lines in the NW and SE cones are about a 
factor of 10 or more, fainter than those in the nuclear re- 
gion using Chandra data. In the NW cone we have found 
the O VII triplet, O VIII RRC, FeXX 3d2p, NeX Lya, 
Mg XI triplet and Si XIII triplet. A similar result is 
found in the SE cone with the detection of C V H7, N 

VII Lya, O VII triplet, O VIII RRC, Ne IX triplet, Fc 
XX, Ne IX ReS, and Mg XI triplet. All of these lines were 
detected in the ACIS / Chandra nuclear region, except O 

VIII RRC. 

At first glance the spectra of the cones look rather dif- 
ferent (see Figure [7]) , the NW cone spectrum shows an 
enhancement around 14A(0.89 keV), which is not present 
in that of the SE cone. To illustrate this effect we have 
constructed two images centered at 0.5-0.7 keV and 0.8- 
1.0 keV, respectively (see Figured]). These spectral inter- 
vals are dominated by the emission of the O VII triplet 
at 0.569 keV and O VIII RRC at 0.871 keV, respectively 
(see Tab. The O VII triplet image shows a mor- 

phology similar to that found in the soft band (< 2 keV, 
see Figure [2]). However, the O VIII RRC image shows 
a bright region to the NW of the nucleus. This can also 




20 30 
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Fig. 7. — Spectral fits ( top panels) and residuals ( bottom panels) 
using low resolution ACIS/ Chandra spectra of the three regions: 
Nucleus (top), cone NW (middle) and cone SE (bottom). 



be noticed in Figure [9] where the radial profiles of the 
O VII triplet and O VIII RRC images are plotted. The 
NW and SE profiles appear consistent in the case of the 
O VII triplet, whereas there is a clear bump at 2.5" in 
the NW side of the O VIII RRC profile compared to the 
SE one. Despite the limited signal-to-noise ratio of the 
spectra from the cones, it seems clear that the O VIII 
RRC is more important in the NW cone than in the SE 
cone (see Figs. Handle. 

5. ORIGIN OF THE SOFT X-RAY EMISSION 

As already said, the soft X-ray spectrum of Mrk573 
is dominated by line emission, as evident from the 
RGS/ XMM-Newton data (Figure gj). Emission lines 
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TABLE 3 

Measured fluxes for EPIC / XMM-Newton and ACIS / Chandra spectra 
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a Units are phem 2 s 1. b This value was obtained by setting the nearby line FeXVII to zero. c These values were computed after 

including NelX HC7 and NeX Ly/3, instead of NelX He<5. 



from H-likc and He-like C, N, O, and Ne and Fe L- 
shell emission line Fe XVII dominate the spectrum, 
which also includes strong narrow RRC O VII and O 
VIII lines. The RRC features of these highly ion- 
ized species are produced when electrons recombinc 
directly to the ground state. They are broad fea- 
tures for hot, collisionally ionized plasma, but are nar- 
row, prominent features for p hotoionized plasma aris- 
ing in low-tem pe rature material (fLiedahl fe Paerelsl 19961 : 
iLiedahll I1999P) . iGuainazzi fe Bianchil ()2007f ) measured 
in the Mrk 573 RGS spectrum the width of the RRC 
lines O VII and O VIII as cr(0 VII) = 3.5±^ eV and 
a(0 VIII) = 4.8 ±3.1 eV, respectively. These values in- 
dicates a relatively cool photoionized plasma. 

Relative emission strength ratios of the 
He-like ls2p x Pi — » Is 2 1 So (r, resonance), 
ls2p 3 P 2 -> Is 2 %, ls2p 3 Pi -> Is 2 % (i, blended 
inter-combinations) and ls2s 3 Si — > Is 2 1 So (f, forbid- 
den) transitions can discriminate between photoion- 
ization, collisional excita ti on, or hybrid environment 
(iPorouet fe Dubaul 120001 : iBautista fe Kallmanl 120001 : 
iKahn et al.l 20021 ). A weak resonance (r) line, com- 
pared to the forbidden (/) and intcrcombination (i) 
lines, corresponds to a pure photoionized plasma. A 
commonly used line ratio is defined as G = (/ + i)/r. 
For Mrk 573 we measured a value of G = 1.4(0.3 - 12) 
for the ion O VII (see Tab. [5] and Figure 3J . This value 
is lower than that expected for a pure photoionized 
plasma (G > 4), but higher than that of collisional 
ionization. Such a value of G would indicate a hybrid 
plasma where collisiona l processes are not negligible 
(jPorquet fe Dubaull2000i although the error bars do not 
allow a secure result. We could only obtain upper limits 
in the intensity of two of the lines needed to compute 
the G ratio corresponding to the other ions (Ne IX and 



N VI) with detected He-like triplet transitions. The 
lower limit of G obtained using Ne IX ((f + i)/r > 2.7) 
is not conclusive although also points to an hybrib 
plasma. However, the use of G ratios has been ques- 
tioned by several autho rs (e.g. IKinkhabwala ct al J [20021 : 
IPorter fe Ferland I2007D . An important contribution of 
photoexcitation would rise the intensity of the resonance 
transition, which yields a decrease of the G ratio relative 
to the pure photoionization case. Kinkhabwala ct al. 
(|2002D proposed the use of the ratio of higher order 
transitions to the forbidden triplet component in Hc-likc 
ions as a good discriminan t between photoex ci tation 
and collisional ionization. [Kinkhabwa la et all (|2002[ ) 
predicted a value of F(0 VII He 7 )/F(0 VII (Hea(f)) 
of 0.017 for collisional ionization. Thus, an excess 
of He7 line relative to the forbidden transition, as 
measured in Mrk 573, F(0 VII He7)/F(0 VII (Hea(f)) 
= 0.31(0.03-1.4) (see Tab. EJ), indicates the impor- 
tance of photoexcitation. 

The intensities of the emission lines Fe XVII L (3d- 
2p) and Fe XVII K (3s-2p) are comparable in the 
case of Mrk 573 (sec Tab. [2]), which could be ex- 
plained as an intermediate temperature plasma (10 eV < 
kT < 500 cV), or under hybrid conditions combining 
collisional and phot oionization equilibrium condit ions 
(jLiedahl et all [1990T) . Similarly, iSako et al.l pOOOT ) ar- 
gued that the dominance of Fe XVII L lines can only 
be explained if photoexcitation by the nuclear radiation 
plays an important role, consistent with our suggestion 
throughout the G ratio. 

From above, we conclude that the emission line diag- 
nostics seems to indicate that the soft-X ray spectrum 
of Mrk 573 as seen by RGS /XMM-Newton data can be 
produced by a plasma dominated by photoionization and 
photoexcitation. Nevertheless, we cannot rule out some 
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Fig. 8.— Images centred at O VII triplet (top) and O VIII RRC 
(bottom) lines. [O III] contours are overlaid. 



contribution of collisional excitation, based primarily on 
the presence of the Fe XVII 3d-2p transition. Unfortu- 
nately, the limited count level of the spectra does not 
permit us to give secur e diagnostics . Deta iled optical 
spectroscopic studies by iFerruit et all (|1999t) found that 
ionizing photons originating in the central source are not 
sufficient to explain the emission line luminosity. They 
suggest that fast shocks, associated with the jet/gas in- 
teractions, might contribute to the gas ionization. 

The ACIS / Chandra spectrum of the nuclear region (1 
arcsec) can be reproduced using the template obtained 
from the RGS/ XMM-Newton analysis. The resulting 
fluxes are consistent with each other. Therefore, it seems 
that most of the emission line fluxes measured in the RGS 
spectra come from the inner central 2 arcsec region and is 
consistent with the photoionization and photoexcitation 
dominated plasma scenario. 

NW cone emission differs from that of the SW cone 
and the nucleus. The O VIII RRC/O VII triplet ratio 
found in the NW cone (F(OVIII RRC)/F(0 VII(r,i,f)) 




2 3 4 5 6 

Radii farcsecs) 

Fig. 9. — Radial profile of the images centred at O VII triplet 
(top) and O VIII RRC (bottom) lines. Blue-dotted line and blue 
dots correspond to the NW cone while black-continous line and 
black dots correspond to the SE cone. 



= O.60 }) i s nominally higher than in the SE cone 
(0.13g;g|) and in the nucleus (< 0.19), although consis- 
tent within the statistica l uncertainties. According to 
iKinkhabwala et ail (|2002f) radiative decay following pho- 
toexcitation dominates the Seyfert 2 spectrum at low 
column densities, whereas recombination following pho- 
toionization dominates at high column densities. Follow- 
ing their prescription, the RRC intensity will increase 
compared to He-like triplet line when column density in- 
creases, which would imply a value of Nh larger for NW 
than for SE. Alternatively, the RRC line could be pro- 
duced by hot collisional ly ionized plasma, although the 
line widths reported by iGuainazzi fc Bianchil (|2007| ) do 
not support this scenario. It is very unlikely that a broad 
feature contributes only in a few percentage of the to- 
tal flux of source. Moreover, the radio maps of Mrk 573 
are not sensitive enough to show great detail concern- 
ing radio emission although two faint blobs are observed 
in the NW cone (jFalcke et al.lll998D . The low signal- 
to-noise ratio of the spectra from the cones prevents us 
from a more detailed study of the line transitions in the 
extended emission. 

5.1. Comparison with photoionization models 

Using ve rsion c08.00 of the Cloudy package (last de- 
scribed by iFerland et al.lll998|) . we attempted to repro- 
duce the observed spectra of the nuclear region and NE 
and SW cones seen in ACIS / Chandra data. 

In these Cloudy simulations we assumed the source of 
ionization to emit as a typical AGN continuum (we used 
the model AGN available in Cloudy) defined by a "big 
bump" of temperature T = 10 6 K, an X-ray to UV ratio 
ctox = —1-15, plus a X-ray power-law of spectral index 
of a = —1.0. A plane-parallel geometry is assumed, with 
the slab depth controlled by the hydrogen column density 
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Fig. 10. — (Top): ACIS / Chandra spectral fit and residuals to 
BLR+NLR Cloudy model of the nuclear region. (Bottom): Same 
for the cone-like regions. Red (black) crosses and red (black) con- 
tinuous line are the NW (SE) cone spectrum and fit, respectively. 



parameter (Nh)- 

Two grids of parameters were constructed to simulate 
the expected BLR and NLR conditions. For each of 
them, a grid of models was simulated by varying the ion- 
ization parameter (U), the density of the material (njj), 
and the hydrogen column density (Nh). U ranges from 
10~ 3 to 10 3 , Nh ranges from 10 20 cm -2 to 10 24 cm" 2 and 
nn from 10 9 cm -3 to 10 11 cm -3 for BLR conditions, and 
from 10 2 cm -3 to 10 4 cm" 3 for NLR conditions. 

In both BLR or NLR models the outputs of each 
Cloudy simulation include a transmitted and a reflected 
emission line spectra. Under the Cloudy terminology "re- 
flected" spectrum refers to the emission escaping into the 
2tt sr subtended by the illuminated face towards the ion- 
izing source and by "transmitted" the emission escaping 
in the opposite direction. Each transmitted or reflected 
simulated emission line spectrum is importe d as additive 
tables ("atable") in XSPEC following iPorter et all 
(|2006l) procedure. All components are also absorbed 
through the Galactic value by using "wabs" on XSPEC. 
In practice, any additive table is included in XSPEC 




3 20 22 
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Fig. 11. — (top): RGS/ XMM- Newton spectra fitted to the best 
fit fou nd for ACIS / Chandra spectrum of the nuclear region (see 
Figure Hot , (bottom): Same than top figure but with the inclusion 
of the thermal model APEC with kT=0.4±0.1 keV. 



as follows: wabs[NH(Gal)](atable{Cloudytable}), 
where "Cloudytable" can be transmitted/reflected for 
BLR/NLR conditions. 

For the nuclear region, neither BLR nor NLR mod- 
els produced any acceptable fit (x 2 = 1.9 and Xr = 1-8, 
respectively). Qualitatively, the issue is that BLR and 
NLR models alone failed in reproducing the relative 
strength between 0.55 keV (the O VII triplet) and 0.9 
keV (the O VIII RRC and/or Ne IX triplet). We have 
found that the result of our Cloudy simulations do not 
change with the value of the volume density and col- 
umn density in the range considered here. The models 
are mostly sensitive to the value of the ionization pa- 
rameter. Besides the two simulations for BLR and NLR 
conditions (hereinafter called BLR and NLR models), we 
have also tried to model the spectrum of the nucleus as 
a combination of both BLR and NLR conditions (here- 
inafter BLR+NLR model). In the BLR+NLR model we 
use them all together, being then two transmitted and 
two reflected emission line spectra. The combination of 
BLR+NLR model produced a good fit (see Figure [T0| 
with Xr = 1-3- The parameters of the best fit are given 
in Table SJ The modelled nuclear spectrum is equally 
dominated by the BLR and NLR emission (BLR is ~54% 
of the emission). We emphasize that in our simulation 
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TABLE 4 

Best fit parameters to a Cloudy BLR+NLR model. Nucleus and Cones. 





Nucleus 


Cones (NW/SE) 




Model param. 


BLR NLR 


NLR 1 


NLR 2 



log(U) 1.23j f| 0.13[| -\l 0.9 ±0.3/0.3 ±0.2 -3*/-3* 

log(n H ) 9.8i°7° 3.14|18 3 * 3. 

log(N H ) [cm 2 ] 20.6|g g 20.6|g | 20* 20* 

Flux(0.5-2.0 keV) (trans.) 35.8 ± 9.6 0.0 (< 6.6) 0.0 (< 5.3) / 9.3 (< 11.6) 3.4 (< 10.9) / 3.1 (< 9.3) 

Flux(0.5-2.0 keV) (reflcc.) 61.1 ± 9.9 92.6 ± 3.7 14.0 (4.6 - 23.0) / 0.0 (< 14.5) 3.8 (< 10.9) / 3.0 (< 9.3) 



Note. — NLR1 and NLR2 refer to high and low ionization phases for the cones using ACIS/ Chandra data. Fluxes in units of 
10 _15 erg cm -2 s _1 . 
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Fig. 12. — 2-D iso-chi-squared flux contours of the normalization 
parameters of the best-fit Cloudy model of the reflected compo- 
nents of the NW (top) and SE (bottom) cones. 



the use of a BLR component does not necessarily implies 
high density material, it rather accounts for ionization 
conditions higher than those of the NLR component. 

In principle we could also expect some contribution 
from collisionaly ionised plasma (see Section [SJ . How- 
ever, we know that this contribution should be small as 
indicated by the G ratios, and discussed in previous sec- 
tion. As an additional check, we used then this best-fit 



solution in order to reproduce the RGS / XMM-Newton 
spectra. We froze all the parameters adding a constant to 
the fit since we were mainly interested to check whether 
this model could reproduce the high resolution RGS spec- 
tra. We preferred this approach of fitting first the low 
resolution spectrum and taking this as a starting point to 
model the high resolution spectra. The final fit is given 
in Figure [TT] (top) with a constant value of 1.2 ± 0.2. It 
gives a good representation of the data although it fails 
in reproducing the FeXVII emission lines. Actually the 
inclusion of a thermal model (APEC) with a kT=0.4±0.2 
keV reproduces better these features (AC ~ 14, see Fig- 
ure [TTJ bottom) . In the later case, it becomes very com- 
plex to discriminate between fits due to the low count 
level present in the high resolution spectrum. Based on 
this result, we went back to the low resolution data and 
try to fit the Chandra nuclear region adding the APEC 
component. Although it does not give formally a better 
fit, we obtained a fraction of the thermal component of 
~ 6% of the total nuclear flux. This results confirms that 
the thermally ionized component is present, although its 
contribution is small. 

We have also modelled the extended emission using 
the grid of Cloudy simulations. In this case we have to 
include the contamination from the PSF wings of the nu- 
clear source. For this, we have computed the expected 
fraction of the nuclear flux contributing to the cone- 
like regions as F n uc2cono =FluxpsF(cone)/ FIuxpsf 
(nucleus). Assuming a Gaussian profile for the Chan- 
dra PSF with a FWHM of 1.2 arcsecs (see Fig- 
ure [TJ), the fraction of nuclear spectrum contribut- 
ing to the cone is -F n uc2conc = 0.006. This factor 
is equivalent to an integrated flux of F(0. 5-2.0 keV)= 
5.06 x 10~ 16 erg cm _2 s _1 , which corresponds to 2.5 and 
3.5% of the NW and SE cone fluxes, respectively. In 
order to account for this nuclear contribution we have 
used the BLR+NLR best fit model found for the nu- 
cleus with all the parameters frozen, and scaled by the 
factor F nuC 2conc- Initially we have tried to fit the cone 
emission by a single phase medium, using Cloudy mod- 
els with NLR conditions and both reflection and trans- 
mission components. Both NW and SE cones have been 
fitted simultaneously, although the normalizations of the 
reflection and transmission components for the two cones 
are let to vary independently The hydrogen column 
density and density of the material have been frozen 
at log(Nn)=20 and log(nn)=3 for simplicity. We recall 
that the simulated spectra by our Cloudy models are 
not sensitive to the value of these parameters within the 
explored range. The best fit show ionization parame- 
ters of log(U)=0.71±0.07 and log(U)=0.13±0.09 for NW 
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and SE cones, respectively. Nevertheless, this fit failed 
to simultaneously reproduce the two spectra, being the 
statistics rather poor (Xt =1-6). A much better fit is 
obtained (x? = 1.1) if we add to this phase (hereinafter 
NLR1) a second one (hereinafter NLR2) with NLR con- 
ditions but a lower value of the ionization parameter 
U. We note that due to convergence problems within 
XSPEC we had to freeze the ionization parameter of the 
second phase. We checked that the best fit is obtained 
for a value of log(U)~ —3, in both regions. The ioniza- 
tion parameters for the NLR1 phases are log(U)=0.9±0.2 
and log(U)=0.3±0.2 for NW and SE cones, respectively. 
These values are consistent with the parameters of the 
previous model. The final fit can be seen in Figure [10] 
and the best fit parameters and fluxes for each model 
component are given in Table 2] In order to show the 
confidence level of the fluxes, Fig. [12] includes the iso- 
chi-squared flux contours of the normalizations of the 
two phases (i.e. NLR1 and NLR2) of the reflected com- 
ponents for the NW (top) and SE (bottom) cones. Sim- 
ilar result is obtained for the transmitted components. 
The NLR1 phase mostly contributes to the X-ray spec- 
trum at the region between 0.8-0.9 keV and 0.53-0.7 keV, 
whereas the NLR2 phase contributes at 0.52 and 0.7-0.85 
keV. The lowest value of the ionization parameter is very 
similar to that required to fit the optical spectrum of 
Mrk 573 (jKraemer et all 120091 ) . In that work, the au- 
thors claimed a three phase component to explain the 
optical emission line spectrum. The low-ionization gas 
accounts for the [OII]AA3727A and [NII]AA6548, 6584A 
emission, whereas the moderately ionized phase accounts 
for the [OIII]AA5007A. Nonetheless, our NLR1 exhibits 
a value of U higher than their highly ionized phase. In 
addition, we expect a contribution of collisionaly ionized 
plasma (see Section[5]) to be present in the extended emis- 
sion. However, given the fact that its contribution to the 
nuclear spectrum is about 6%, we assumed its contribu- 
tion to the extended emission to be equal or smaller than 
in the nuclear case and did not try any fit given the low 
count level in these regions. 

In terms of flux, the NLR1 phase is the 64% of the ex- 
tended emission flux. Its contribution in the NW cone is 
higher than in the SE cone. Moreover, the reflection com- 
ponent dominates the NW cone whereas the transmission 
component dominates the SE one. Tentatively, one could 
attribute this result to an orientation effect, being the 
NW cone located in the farthest side and the opposite 
for the SE one. In deed, 2-D spectroscopic observations 
(|Ferruit et allll999T ) indicates that at least part of the NW 
cone ionized gas is red-shifted with respect to the sys- 
temic velocity. This could indicates that NW cone axis 
could be oriented behind the sky plane, although close 
to it. This orientati on apparently contradicts t he ori- 
entation proposed by iTsvetanov fc Walsh] (|1992[ ) based 
on reddening measurements, although we remark that 
their findings refer to the orientation with respect to the 
galaxy disk, not to the sky plane. Remarkably, the NLR2 
phase accounts for the same flux in the NW and SE cones 
and nearly equal relative contribution of the reflection 
and transmission components. Thus, this low ionization 
phase seems to be uniformly distributed along the cone 
area. 

At this point, we may conclude that soft extended X- 
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Fig. 13. — Variation of the ratio of soft X-ray to [O III] brightness 
profile in arbitrary units. The [O III] brightness profile variation is 
also ovcrplotcd as a green histogram-like line. Positive x-axis val- 
ues indicates the NW orientation and negative values the SE one. 
Horizontal dashed lines represent predictions from simple photoion- 
ization models at various values of the ionization parameter. 

ray emission is then powered by features coming from 
ionized gas in two phases under NLR-like conditions. 
Furthermore, one may question: is the [O III] extended 
emission structure powered by the same mechanism? 
The similarity between the [O III] structures and the 
soft X-ray emission points to a common origin for both 
components. We have investigated the radial variation of 
the ratio between the soft X-ray and the [O III] line emis- 
sion. We represented in Figure Q2] the variation of the 
brightness ratio along the axis of the cone (PA = 122°). 
The brightness profiles were extracted using the IRAF 
pvector task. Before obtaining the ratio we convolved 
the [O III] profile to obtain the resolution of X-ray data. 
As can be seen from Figure [T3] the ratio soft-X/[0 III] 
presents a non-uniform variation showing a maximum at 
the nucleus; it then drops dramatically at the position 
of the [O III] arcs and returns at roughly half of the nu- 
clear value outside the [O III] arcs. A similar behaviour, 
na mely a small y a riatiq n of the ratio, has been reported 
bv lBianchi et~aTl (120061) for the case of NGC 3393. We 
have also compared the radial variation with predictions 
from photoionization models in Figure [T3l We have used 
simple models assuming single plane-parallel slabs, con- 
stant density and radiation bounded clouds. The soft 
X-ray emission has been taken as the sum of the pre- 
dicted values for the most intense features identified in 
our X-ray spectra. We have scaled the model predic- 
tions to the observed nuclear values from the model with 
logU = 0.5, which is close to the best fitting value de- 
rived above. The predictions for different values of the 
ionization parameters are represented in Figure 1131 It 
can be seen that a variation of Ut by one order of mag- 
nitude is needed to reproduce the observed variation in 
the brightness ratio from the nucleus to the arcs, which 
could be attributed to a combination of radiation dilution 
plus density enhancements at the arc positions. Note, 
that these variations are qualitatively in agreement with 
our Cloudy model simulations in which two different U 
values are needed for SE and NW cones. A baseline 
m odel where the densit y decreases as r~ 2 , as proposed 
by iBianchi et al.l (]2006l ). is compatible with our results, 
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Fig. 14. — High spatial resolution nuclear SED of Mrk 573 fitted 
with the clumpy torus models. Solid and dashed lines are the best 
fitting model and that computed with the median of each of the 
six par ameters that describe the models (sec Ramos Almeida ct al. 
1200911 ). The shaded region indicates the range of models compat- 
ible with a 68% confidence interval around the median. 

except for the regions close to [OIII] arcs. It is very 
likely that simple photoionization models are not ade- 
quate to reproduce the observed ionization variations, 
although a more sophisticated treatment must wait until 
higher quality X-ray spectroscopic observations become 
available. The fact that ionization- and matter-bounded 
clouds arc likely constituents of the NLR has not been 
explored to explain emission lines in the soft-X range. 

6. THE NUCLEAR SPECTRAL ENERGY DISTRIBUTION 

Another important feature of the unified model is the 
optically thick torus. Under this scheme, Type-2 Seyferts 
like Mrk 573 are obscured due to this material located 
along our line of sight. The best way to study and char- 
acterize the molecular torus is by modeling the nuclear 
spectral energy distribution (SED) of the sources. The 
near- and mid-IR nuclear emission of Scyfert galaxies is 
attributed to the reprocessing of the UV/X-ray nuclear 
radiation by the toroidal dusty structure. For this reason, 
the infrared range is key to put constraints on torus mod- 
eling. However, in comparing the predictions of any torus 
model with observations, the small-scale torus emission 
must be isolated, in order to avoid contamination from 
the host galaxy. For this reason, it is important to use 
high angular resolution data when trying to model the 
torus emission. 

We have tried here to explain the optical/infrared nu- 
clear SED of Mrk 573 constructed with high spatial res- 
olution data (see Tabic [T]) using an interpolated ver- 
sion of t he recent mod e ls for the clumpy torus sce- 
nario by iNenkova et~aT1 (|2008alfbh . We have searched 
for the best fitting models using the Bayesian inference 
tool B ayesClumpy (jAsensio Ramos fe Ramos Almeidal 
120091) The results of the fit are shown in Figure HU 

Indeed, the nuclear SED of Mrk 5 7 3 has b een previ- 
ously fitted bv IRamos" Almeida et all (|2009bt l using this 
set of models and tools, although the data points be- 
low 1 [im were not included in their analysis. For using 
the optical photometry derived in this work, the clumpy 
torus model fit needs an additional extinction factor, 
which is included as a foreground extinction. The derived 
median value is Ay = 7.3 ±0.5 mag, which can be trans- 



lated to a column density of Nh = 1.39 x 10 22 cm" 2 . 
This value is nicely consistent with that derived from 
the optical colour-colour maps (see Sect. [3]). 

The results of the fitting process are the probability 
distributions for the free parameters that describ e the 
clumpy models (see IRamos Almeida et all 1200 9b). For 
Mrk 573, the median values of the parameters corre- 
spond to a torus width of 30°, N = 4 ± 1 clouds in 
the equatorial direction, optical depth per single cloud 
t v ~ 60, and A^ os - 300 mag (equivalent to N[j os ~ 
5.7 x 10 23 cm" 2 ). These torus parameters are similar to 
th ose derived without including the optical data points 
bv IRamos Almeida et all (|2009b[ ). 

We have determined the torus luminosity integrating 
the corresponding emission from the torus model corre- 
sponding to the median value of the probability distribu- 
tion of each parameter (dashed line in Figure [T4")) . The 
resulting value is L£°; = 7.4 x 10 43 ergs -1 . The clumpy 
model fit yields the bolometric luminosity of the intrinsic 
AGN, L^f N = 5.2 x 10 44 ergs" 1 (the bolometric lumi- 
nosities are good to a factor of 2). Combining this value 
with the torus luminosity, we derive the reprocessing effi- 
ciency of the torus, which for Mrk 573 results to be quite 
low, about 14%. 

Moreover, we have derived the X-ray luminosity from 
the hard X-ray part, assuming a power law with pho- 
ton index 1.8. The range from 2 to 6 keV was used 
to avoid, on the one hand the soft range, which is at- 
tributed to emission lines, and the Fc K-line on the 
other. The galactic absorption plus an intrinsic absorp- 
tion were included in the fitting process. The intrinsic 
absorption results in a value of Nh = 7lg x 10 22 cm" 2 , 
which is compatible with the foreground extinction in- 
ferred from the colour maps and the clumpy torus mod- 
eling. Thus, the absorption corrected X-ray luminosity 
is L 2 -io kcv = 1.7 x 10 ergs -1 . This value is quite low 
when compared to the optical/infrared luminosity repro- 
cessed by the blocking torus. However, it can be recon- 
ciled by taking into account that Mrk 573 has been clas - 
sified as a Compton-thick AGN (jGuainazzi eTaDEool), 
and its X-ray luminosity has to be corrected by a large 
factor to obtain the intrinsic luminosity, since the column 
density is Nh > 1.6 x 10 24 cm -2 (such high value of the 
optical depth is also consistent with the value derived 
above from the clumpy torus modeling). iPanessa et al.l 
(|2006h derived a factor of 60 comparing a small sam- 
ple of Compton-thick Typ e-2 Seyferts with a sample of 
Type-1 Seyferts, whereas iCappi et al.1 (|2006| ) derived a 
factor abo ut 100. These values are in contrast with 
the work of lGon zalcz- M"artin et al.l (|2009f ). who obtained 
a value of 42 using a sample of LINERs. In addi- 
tion, we have to transform from X-ra y to bolometric 
lumin osity mu ltiplying by a facto r 30 (jRisaliti fc Elvisl 
120041 : see also IPanessa et all 120061 ) . Thus, we obtain a 
value for the bolometric AGN luminosity in the range 
Lagn = 3.1 — 5.1 x 10 44 ergs -1 , which is in nice agree- 
ment with the value derived from the t orus reprocess- 
ing, g iven the uncertainties involved. iKraemer et al.1 
(|2009f l derived a higher value for the bolometric lumi- 
nosity (3.2 x 10 45 ergs' 1 ) based on the [OIV] 25.9^m 
luminosit y, as measured by the Spitzcr/IRS spectrum of 
Mrk 573 (jMelendez et al.ll2008al lbh. In any case Mrk 573 
seems to be radiating near to the Eddington limit as- 
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suming the black hole mass to be around 2 x 10 7 M Q 
(|Bian fc Gull2007f) . This result adds support to the re- 
cl assification of Mrk 573 as a h idden narrow-line Seyfert 

1 (jRamos Almeida et alJl2008f >. 

7. CONCLUSIONS AND OVERALL PICTURE 

Mrk 573 is a nearby optically classified Type-2 Seyfert, 
well-known for its extended circumnuclear emission-line 
regions. It is this extension and the proximity of 
the source that convert Mrk 573 as one of the ideal 
cases to study this emission commonly found in Type- 

2 Seyfert galaxies. We combine RGS/ XMM-Newton and 
ACIS / Chandra to achieve high spectral and spatial reso- 
lution in order to disentangle the emission mechanism of 
this extended emission. We also used optical and near- 
IR HST data in order to compare with the X-ray data. 
The main results are: 

• The soft X-ray emission is very complex, resem- 
bling that of the [O III] emission, as already re- 
ported. What constitutes a new result is that espe- 
cially the NW structure is also very similar to that 
of [O III] /Ha emission. This suggests the same 
origin for the emission lines at optical and the soft 
X-ray ranges. 

• Through X-ray spectroscopic analysis wc have 
found that plasma excitation mechanism in the nu- 
clear spectrum is mainly driven by photoionization 
from the central source, including a strong contri- 
bution from photoexcitation. A small contribution 
of collisionally ionized plasma is also needed to ex- 
plain the emission line ratios shown by RGS spec- 
tra. This conclusion also agrees with the proposed 
Cloudy simulations since the spectra could be inter- 
preted as the combination of two different phases 
of Cloudy models, with two different ionization pa- 
rameters log(U)=1.23 and log(U)=0.13. 

• Based on the ACIS/ Chandra images and radial 
profiles along O VII triplet and O VIII RRC for 
cone-like structures we showed that O VIII RRC 
could be more relevant toward the NW cone, al- 
though the line ratios are formally compatible after 
including the error bars. 

• We have successfully modelled the cone-like emis- 
sion using Cloudy simulations corresponding to 
two phases of NLR conditions. The first phase 
shows different values of ionization parameter and 
different contributions of the reflected and trans- 
mitted components for the NW (log(U)=0.9, re- 
flected dominated) and SE ( log(U)=0.3, transmit- 
ted dominated) cones. The second is an homo- 
geneous phase with a lower ionization parameter 
(log(U)=-3) and the same contribution of reflected 
and transmitted componcts for NW and SE cones. 

• We have found a good agreement between the AGN 
bolometric luminosity derived from the hard X-ray 



luminosity and the one derived from the modeling 
of the optical/infrared SED with a clumpy torus 
model. 

• From the extinction maps we have found that a 
dust lane crosses the nucleus in the N-S direc- 
tion. This could be, in projection, perpendicular 
to the direction of the cone-like structure. The 
amount of extinction we have derived from this 
map is also consistent with that derived from the 
SED modeling. However, this must be related with 
the extended extinction since, in the inner parts, 
the AGN is hidden by a column density on the 
Compton-thick regime (Nh > 1-6 x 10 24 cm — 2). 



Note added in manuscript: After this paper was 
submitted to the jo urnal a paper was published by 
iBianchi et "all (|2010t l with data in common with our 
work. Most of their results are consistent with ours, al- 
though different approaches were used. They fitted the 
XMM-Newton /RGS spectra using Cloudy photoioniza- 
tion models. Their best fit was obtained by an hybrid 
model -photoionization + collisional excitation- where 
the collisional phase contributes 1/3 of the flux in the 
band 0.5-0.8 keV, consistent with our work. Moreover, 
they claimed the need of two photoionization phases (log 
U = 0.3 and 1.8) to explain the ACIS/ Chandra spectrum 
in the range 0.4-7 keV, which is also in agreement with 
our results. The results on the extended emission cannot 
be compared since they do not distinguish between the 
two cone-like structures, extracting the spectrum from a 
circumnuclear annulus. 
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